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1.  INTRODUCTION 


There  are  a  number  of  munition-related  problems  of  interest  to  the  defense  community  in  which  the 
presence  of  energetic  materials  is  integral.  One  such  area  of  investigation  is  that  of  explosive  ordnance 
disposal  (EOD)  render-safe  procedures,  in  which  a  potentially  armed  and  thus  dangerous  fuze  is  disabled 
in  situ,  in  order  to  render  it  safe  for  subsequent  handling  and  removal.  One  class  of  render-safe  methods 
might  be  termed  "mechanical,"  in  that  it  relies  on  a  measured  use  of  mechanical  force  to  sever  or  disrupt 
the  explosive  initiation  train  within  the  fuze. 

One  such  mechanical  technique  to  render  a  fuze  safe  is  to  impact  it  with  a  hardened  projectile  at  an 
oblique  angle,  for  example,  perpendicular  to  the  fuze's  axis  of  symmetry.  With  the  proper  selection  of 
impact  location  and  speed,  the  mechanical  damage  induced  by  the  impact  can  be  made  sufficiently 
brisant  to  isolate  the  high-explosive  primer  material  from  subsequent  stages  of  the  high-explosive  train. 
However,  the  impact  should  not  shock  the  fuze  to  an  extent  which  causes  the  explosive  material  to 
initiate.  The  EOD  community  currently  employs  this  mechanical  method  in  certain  render-safe 
procedures. 

The  physical  process  of  rendering  a  fuze  safe  in  situ  can  be  quite  hazardous.  In  addition,  certain 
types  of  foreign-made  ordnance  are  of  limited  supply,  and  not  available  for  parametric  testing.  For  these 
reasons,  it  is  recognized  that  numerical  simulation  may  assist  in  the  development  of  render-safe 
procedures.  To  this  end,  this  report  discusses  an  investigation  into  the  suitability  of  employing  the 
EPIC92  hydrocode  in  the  modeling  of  an  oblique  fuze  impact,  a  mechanical  render-safe  procedure.  The 
corresponding  experiment  was  performed  after  the  simulation  work  was  completed,  which  provides  a 
direct  comparison  for  the  simulation  results.  Another  advantage  of  conducting  the  impact  experiment 
after  the  numerical  analyses  is  in  determining  the  best  time  to  radiograph  the  fuze.  Upon  introducing  the 
specifics  of  the  idealized  problem  to  be  studied,  the  problem  geometry,  constitutive  properties,  and 
boundary  conditions  of  the  problem  are  described.  The  results  of  finite  element  numerical  simulations 
and  experiment  are  presented  and  compared.  Finally,  an  elastic-plastic  analytical  solution  is  developed, 
which  sheds  insight  into  fuze  render-safe  mechanics,  and  is  compared  with  the  hydrocode  results. 
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2.  PROBLEM  DESCRIPTION 


Current  work  focuses  on  investigating  the  "structural"  damage  (i.e.,  the  mechanical  response) 
imparted  to  an  inert  fuze  simulant  as  a  result  of  an  oblique  impact,  rather  than  establishing  the  impact 
conditions  that  are  necessary  to  render  safe  an  actual  fuze  device.  Only  after  a  proper  understanding  of 
the  mechanical  response  is  obtained  may  future  work  address  the  complete  fuze  render-safe  problem  by 
replacing  the  inert  components  of  the  fuze  simulant  with  an  explosive  material.  Hereafter,  the  fuze 
simulant  will  simply  be  referred  to  as  the  fuze. 

A  typical  fuze  geometry  consists  of  a  right-circular  cone  that  threads  into  the  ordnance  at  its  base 
(Figure  1).  The  cutaway  view  also  illustrates  some  of  the  important  internal  fuze  components  such  as 
the  fuing  pin,  detonator,  and  booster.  Although  these  components  are  critical  for  proper  operation  of  the 
fuze  under  normal  impact  conditions,  replacing  the  "fine"  internal  structure  of  the  fuze  with  a 
homogeneous  core  material  will  simplify  the  oblique  impact  problem. 

The  fuze  selected  for  analysis  is  comprised  of  an  axially  symmetric  6061-T6  aluminum  sheath 
surrounding  a  softer  brass  core.  The  aluminum  sheath  is  2.40  in  (61.0  mm)  in  diameter  at  the  base,  and 
tapers  bilinearly  from  base  to  apex.  The  fuze  height  is  3.72  in  (94.5  mm).  The  projectile  consists  of  a 
solid,  hard  steel,  right-circular  cylinder,  1  in  (25.4  mm)  in  diameter  and  3  in  (76.2  mm)  long,  with  a  mass 
of  300  g.  It  is  designed  to  strike  the  fuze  normal  to  its  axis  of  symmetry  during  standard  render-safe 
procedures.  The  initial  configuration  as  modeled  is  shown  in  Figure  2,  including  a  typical  hit  location 
and  impact  velocity.  The  specific  engagement  configurations  for  the  three  numerical  simulations  and 
single  experiment  were  provided  by  Patel  and  Gold  (1993)  (see  Sections  3.2  and  5). 

It  should  be  pointed  out  that  in  an  actual  oblique  impact  render-safe  procedure,  the  barrel  of  the 
gun  which  launches  the  projectile  is  kept  in  close  proximity  to  the  target  fiize.  This  actual  configuration 
has  the  effect  of  providing  lateral  confinement  on  the  projectile.  However,  both  experiment  and 
simulations  are  performed  with  a  free-flying  projectile,  because  a  more  proximate  gun  position  might 
obscure  the  photographic  coverage  with  propellant  combustion  products.  This  lack  of  lateral  confine- 
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ment  has  two  effects:  it  permits  projectile  rotation  after  impact,  and  it  reduces  the  accuracy  of  the  aim 
in  the  experiment. 


3.  NUMERICAL  APPROACH 

A  number  of  hydrocodes  (for  example,  CTH  (McGlaun  et  al.  1987),  SPH  (Monaghan  1982),  CALE 
(Tipton  1991),  EPIC  (Johnson  1977;1992a,b),  MESA  (HoUan  et  al.  1991),  HULL  (Matuska  1984),  and 
DYNA  (Hallquist  and  Benson  1986))  are  available  to  numerically  simulate  the  render-safe  impact  event 
in  question.  The  term  "hydrocode"  describes  that  class  of  continuum  mechanics  based  numerical  codes 
capable  of  solving  wave  propagation  problems.  This  investigation  uses  the  Lagrangian-based  explicit 
EP1C92  (Elastic -Plastic  Impact  Calculations)  hydrocode  (Johnson  et  al.  1992a,b)  because  of:  1)  previous 
impact  modeling  experience  with  EPIC92  and  other  similar  Lagrangian-based  hydrocodes  (Segletes 
1 990),  and  2)  EPlC's  capabilities  for  modeling  ignition  and  detonation  phenomena  in  explosive  materials 
(Johnson  et  al.  1992a,b;  Jones  and  Zerilli  1993).  This  explosive  modeling  capability  will  permit  future 
efforts  to  model  a  more  realistic  fuze  configuration  while  at  the  same  time  retain  compatibility  with  the 
current  series  of  simulations. 

Since  t^e  three-dimensional  dynamic  equations  of  motion  for  the  oblique  fuze  impact  problem  are 
intractable  to  solve  in  closed  form,  finite  element  methods  are  employed  to  obtain  approximate  solutions 
to  the  problem.  A  complete  description  of  any  boundary  value  problem  in  continuum  mechanics 
involves  specifying  the  geometry,  boundary  conditions,  constitutive  equations,  and  failure  or  instability 
criterion.  The  fmite  element  solution  for  the  unknown  field  variables  (typically  displacements)  over  the 
domain  of  interest  proceeds  by  using  the  method  of  weighted  residuals  or  through  the  development  of 
a  variational  principle  for  the  problem  (see  Zienkiewicz  (1977)  for  a  good  review  of  finite  element 
methods).  As  described  in  Zienkiewicz  (1977),  a  global  stiffnes  matrix  is  assembled  at  each  time  step 
to  obtain  a  finite  element  solution.  Alternatively,  the  equations  of  motion  are  directly  integrated  in 
EPIC92,  thus  circumventing  the  need  for  the  formation  of  a  stiffness  matrix  (Johnson  1977). 

3.1  Geometry.  A  typical  fuze  geometry  consists  of  a  right-circular  cone  that  threads  into  the 
projectile  (Figure  1).  The  cutaway  view  also  illustrates  the  important  internal  fuze  components  such  as 
the  firing  pin,  detonator,  and  booster.  Although  these  complex  components  are  necessary  for  operation 
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of  the  fuze  under  normal  impact  conditions,  we  replace  them  with  a  homogeneous  material  to  simplify 
the  analysis.  Thus,  the  fuze  model  consists  of  an  axially  symmetric  outerbody  that  encases  an  inner 
homogeneous  core  as  shown  in  Figure  2.  The  base  of  the  fuze  is  2.40  in  (61  mm)  in  diametei;and  the 
fuze  height  is  3.72  in  (94.5  mm)  (Figure  2).  The  projectile  consists  of  a  solid  right-circular  cylinder, 
1.0  in  (25.4  mm)  in  diameter  and  3.0  in  (76.2  mm)  long,  and  impacts  the  fuze  normal  to  its  axis  of 
symmetry  during  EOD  render-safe  procedures. 

3.2  Boundary  Conditions.  The  boundary  conditions  are  usually  of  mixed  type  and  are  given  in 
terms  of  displacements,  forces,  or  their  first  time-derivatives.  The  fuze  is  assumed  to  be  rigidly  fixed  at 
the  base  so  that  all  x-,  y-,  and  z-displacements  vanish  on  the  base  surface  (Figure  2).  We  believe  this 
assumption  will  not  severely  affect  the  results  because  the  base  of  the  actual  fuze  is  tightly  screwed  into 
the  ordnance  (Figure  1).  Since  the  core  material  is  press-fit  into  the  fuze  outerbody,  in  the  simulation 
we  assume  that  the  homogeneous  brass  core  material  is  bonded  to  the  fuze  outerbody  (Figure  2).  The 
projectile  impact  velocity  is  198  m/s,  and  numerical  solutions  are  sought  for  two  impact  points  along 
the  vertical  fuze  axis  (Figure  2).  One  impact  point  is  1.25  in  (31.75  mm)  below  the  top  of  the  fuze,  and 
the  other  impact  point  is  0.50  in  (12.7  mm)  above  the  base  of  the  fuze;  the  bottom  edge  of  the  projectile 
impacts  the  fuze  at  these  positions. 

Finally,  we  assume  that  the  contact  surface  between  the  projectile  and  the  fuze  is  frictionless.  It 
is  difficult  to  estimate  the  severity  of  this  assumption  because  of  the  lack  of  dynamic  friction  data  for  these 
materials.  However,  friction  effects  dominate  the  interactions  of  bodies  in  very  low-speed  or  quasistatic 
contact  problems.  Thus,  in  the  late  stages  of  the  impact  event,  as  the  projectile  decelerates,  friction  will 
play  a  more  important  role  in  the  interactions.  A  three-dimensional  slide  line  algorithm  permits  the 
Lagrangian-based  hydrocode  to  accommodate  the  severe  distortions  encountered  in  impact  problems. 
This  algorithm  typically  involves  defining  a  master  and  slave  surface,  so  as  to  prevent  element 
interpenetration.  A  master  surface  and  slave  surface  exist  for  both  the  projectile  and  the  fuze  in  this 
analysis. 

3.3  Constitutive  Equations.  The  fuze  outerbody  material  consists  of 6061 -T6  aluminum  (EPIC92 
material  model  23)  that  encases  an  inner  core  of  gun-cartridge  brass  (EPIC92  material  model  2)  (Figure 
2).  The  projectile  material  consists  of  4340  steel  (EPIC92  material  model  9).  The  initial  simulation 
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results  reveal  that  the  projectile  deforms  more  than  that  anticipated  (Patel  and  Gold  1993).  In  subsequent 
analyses,  we  use  a  harder  steel  projectile,  with  a  hardness  of  RockweU  C  (Rc)  Rc57.  The  properties  of 
this  material  are  based  upon  a  modification  of  an  Rc50,  S-7  tool  steel  (EPIC92  material  model  10)  with 
a  modified  initial  yield  stress;  the  magnitude  of  the  modified  yield  stress  is  based  upon  an  empirical 
correlation  that  exists  between  Rockwell  hardness  and  yield  strength  in  various  metals  (see  Oberg,  Jones, 
and  Horton  1984).  The  fuze  deformation  mechanisms  observed  for  projectiles  with  these  hardness 
values  (Rc30  and  Rc57)  are  very  different,  and  therefore  will  be  discussed  and  compared  in  more  detail 
in  Section  4. 

The  Johnson-Cook  viscoplastic  constitutive  model  is  used  in  the  analysis  together  with  a  von 
Mises  initial  yield  condition  (Johnson  and  Cook  1985;  Johnson  et  al.  1992a,b).  An  isotropic  hardening 
rule  governs  subsequent  yield  surface  behavior.  The  equivalent  uniaxial  stress,  o ,  is  defined  as 
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in  which  a'.,  is  the  deviatoric  stress  tensor.  The  Johnson-Cook  model  incorporates  strain,  strain  rate, 
and  thermal  softening  effects  where  the  corresponding  parameters  functionally  appear  in  multiplicative 
form  as 
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in  which  e  is  the  equivalent  time-dependent  uniaxial  plastic  strain  defined  by 
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and  dtj’.  =  ttejj  -  .  D,.^,  is  the  fourth-rank  elastic  modulus  tensor. 


is  the  dimensionless  plastic 


strain  rate  for  e„  =  is',  and  T*  is  the  homologous  temperature.  The  material  constants  A,  B,  C,  n,  and 
m  are  determinable  from  material  tests  conducted  at  different  strains,  strain  rates,  and  temperatures 


(Johnson  el  al.  1992a,b).  The  equivalent  adiabatic  flow  stress  versus  equivalent  plastic  strain  plots  at 
strain  rates  of  1,  100  and  10,000  s  '  are  derived  from  the  EPIC92  material  library  constants  (Figure  3). 
All  of  the  materials  (except  for  brass)  exhibit  work-softening  behavior  which  is  a  result  of  adiabatic 
thermal  softening  of  dynamically  strained  material,  rather  than  strength  loss  due  to  microcracking  or 
other  damage  mechanisms. 


3.4  Damage  and  Erosion  Criteria.  The  damage  criterion  in  EPIC92  is  assumed  to  be  a  scalar 
function  of  the  incremental  (cumulative)  equivalent  plastic  strain,  Ai',  normalized  with  respect  to  the 
strain  at  fracture,  (Johnson  et  al.  1992a,b): 


D  = 


(4) 


E 

where  e'  takes  on  a  similar  multiplicative  form  to  that  of  Equation  (2), 


f 

£  = 

Di  +  D2 

1  D4  In  E 

1  -1-D5T* 

L  J 

. 

(5) 


in  which  o’  =  /  o  is  the  pressure-stress  ratio,  is  the  pressure  (mean  stress),  and  through  Dj 

are  independent  material  constants  determined  from  fracture  experiments.  If  D  >  1  in  Equation  (4) 
in  a  given  finite  element,  the  element  can  only  sustain  additional  hydrostatic  compressive  stress  but  not 
shear  or  tensile  stresses.  Unlike  other  finite  element  codes  which  can  model  the  propagation  of  discrete 
fracture  surfaces  in  the  finite  element  mesh  (Ingraffea  1989),  EPIC92  cannot  simulate  such  behavior. 
EPIC92  has  the  capability  of  handling  the  severe  distortions  encountered  in  penetration  problems  using 
an  element  erosion  algorithm.  The  algorithm  eliminates  element  volume,  but  retains  element  mass  when 

a  critical  value  of  i ’’  is  reached.  In  the  current  analysis,  we  consider  that  element  volume  erodes  when 
e”  =  1.5. 


4.  NUMERICAL  RESULTS 


The  three-dimensional  finite  element  mesh  used  to  model  the  fuze  consists  of  544  symmetric 
brick  finite  elements;  each  brick  element  contains  24  tetrahedral  elements  that  are  formed  by  linking  four 
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constant  strain  triangles  (CSTs)  (Figure  4).  There  are  a  total  of  13,056  tetrahedral  elements.  The  low- 
order  displacement  field  within  each  CST  element  necessitates  the  use  of  more  elements  in  regions  of 
high-deformation  gradient.  One  advantage  of  conducting  the  numerical  analyses  prior  to  the  impact 
experiment  is  in  determining  the  best  time  to  radiograph  the  fuze.  The  results  of  the  three-dimensional 
EPIC92  numerical  simulations  appear  in  Figures  5, 6,  and  7.  Figure  5  illustrates  the  "soft"  projectile 
(Rc30)  high-impact  point  results.  Figures  6  and  7  illustrate  the  "hard"  projectile  (Rc57)  high-impact  and 
low-impact  results,  respectively. 

4. 1  Event  1.  "Soft"  Projectile.  Hi eh-Impact  Point.  The  "soft"  projectile  impacts  the  fuze  ogive 
normal  to  the  fuze  axis  of  symmetry  at  a  position  located  1.25  in  (31.75  mm)  below  the  top  (measured 
along  the  fuze  axis)  of  the  fuze.  The  projectile  plastically  deforms  and  rotates  clockwise  as  it 
progressively  slips  towards  the  fuze  apex  (projectile  rotation  sense  is  relative  to  that  shown  in  Figures 
5, 6,  and  7).  The  projectile’s  slipping  motion  scours  a  continuous  groove  of  elliptical  form  into  the  fuze’s 
surface.  However,  the  groove  is  not  so  deep  so  as  to  expose  the  inner  brass  core  of  the  composite  fuze. 
After  approximately  420  ps  of  continuous  contact,  the  projectile  slips  off  the  end  of  the  fuze  and 
continues  rotating  clockwise.  After  1 ,200  |is,  the  rear  edge  of  the  projectile  rotates  sufficiently  to  hit  the 
fuze  in  a  secondary  impact  event;  this  impact  causes  only  minor  additional  fuze  deformation.  Because 
the  projectile  sustains  large  permanent  deformations  in  our  simulations,  while  Patel  and  Gold  (1993) 
have  observed  that  the  projectile  is  relatively  undeformed  after  impact,  subsequent  simulations  use  a 
harder  steel  projectile  on  the  order  of  Rockwell  C  57.  The  "hard"  projectile  material  properties  are 
essentially  that  of  an  S-7  tool  steel  (see  Section  3.3,  Constitutive  Equationsl. 

4.2  Event  2.  "Hard"  Projectile.  High-Impact  Point  Under  identical  impact  conditions,  and  in 
contrast  to  the "  soft"  projectile  impact,  the  "hard  ’’  projectile  impact  is  relatively  slip-inhibited  and  results 
in  a  greater  penetration  depth  into  the  fuze  than  with  the  "soft"  projectile.  The  "hard"  projectile  slips  a 
negligible  distance  in  the  vertical  direction,  however,  and  rotates  counterclockwise  rather  than  clock¬ 
wise,  which  can  be  seen  by  comparing  Figures  5  and  6  at  an  elapsed  time  of  200  ps.  As  in  the  "soft" 
impact  event,  the  "hard"  impact  does  not  expose  the  brass  core  of  the  composite  fuze.  However,  the 
"hard"  projectile  causes  more  deformation  in  the  fuze  than  the  "soft"  projectile;  a  comparison  of 
deformed  meshes  at  300  ps  in  Figures  5  and  6  shows  a  greater  degree  of  bending  in  the  fuze  impacted 
by  the  "hard"  projectile.  Hence,  it  appears  that  a  "soft"  projectile  impact  will  cause  less  fuze  damage 
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due  to  slip  and  permanent  projectile  deformation;  the  "hard”  projectile  causes  more  bending  in  the  fuze 
and  the  projectile  sustains  less  permanent  deformation.  The  contact  duration  for  the  "hard"  projectile 
impact  is  about  500  |is,  which  exceeds  the  durations  of  the  other  two  impact  events. 

Equivalent  plastic  strain,  i'’,  contours  aid  in  visualizing  the  permanent  deformation  history  of 
the  fuze  (Figure  8).  Equivalent  plastic  strain  is  a  scalar  representation  of  the  tensorial  strain  state.  The 

e  =  0.20  contours  develop  asymmetrically  beneath  the  edge  of  the  projectile  after  20  )is  (Figure  8a). 
Plastic  strain  then  spreads  through  the  aluminum  outerbody  and  reaches  the  brass  core  at  40  [xs.  With 
increased  projectile  penetration  and  fuze  bending,  plastic  deformation  also  appears  on  the  back  side  of 

the  fuze,  and  achieves  e'’=0.20  atl20|is.  The  plastic  zones  eventually  merge  to  form  a  plastic  "hinge" 

through  the  cross  section  of  the  fuze.  The  projectile  also  deforms  permanently,  and  i"  =  0.20  is  reached 
beneath  the  impact  edge  after  220  ps.  With  continued  bending  of  the  fuze,  the  plastic  hinge  widens 
somewhat,  but  then  remains  steady  throughout  the  remainder  of  the  impact  process. 

4.3  Event  3.  "Hard"  Projectile.  Low-Impact  Point.  In  contrast  to  the  "hard"  high-impact  event, 
the  "hard"  low  impact  causes  much  less  fuze  bending,  primarily  because  the  low  impact  occurs  near  the 
thick,  rigidly  fixed  fuze  base  where  the  bending  stiffness  is  large  (Figure  7).  Furthermore,  because  the 
impacted  fuze  surface  is  less  oblique,  the  entire  surface  of  the  projectile  contacts  the  fuze  after  20  ps.  This 
reduces  both  projectile  rotation  and  slip  along  the  fuze  surface.  The  equivalent  plastic  strain  contours 

are  nearly  symmetric  about  the  contact  region  (Figure  8b).  In  contrast  to  the  high-impact  event,  the  e 
=  0.20  contours  do  not  extend  to  the  back  side  of  the  fuze,  since  there  is  not  as  much  fuze  bending.  There 
is  also  a  reduction  in  the  contact  time,  and  after  280  ps  the  projectile  bounces  off  the  fuze  in  elastic 
rebound.  The  impact  simulation  results  are  summarized  in  Table  1. 

Table  1.  Summary  of  Impact  Simulation  Results. 
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"•oft,"  high 
impact 

Event  2. 
"hard,"  high 
impact 
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"hard,"  low 
impact 
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5.  IMPACT  EXPERIMENTS 


In  the  experiment  a  300-g  cylindrical  steel  projectile  is  launched  at  208  m/s  from  a  37-mm  gun 
with  a  barrel  diameter  of  1.090  in  (27.7  mm).  A  schematic  showing  the  relative  positions  of  the  fuze, 
projectile  breakwire  fixture  for  velocity  measurement,  x-ray  tube  and  film  plate,  and  16-mm  Hycam 
camera  is  given  in  Figure  9.  The  propelling  charge  consisted  of  28  g  of  M2-165-mm  gun  propellant  and 
was  ignited  using  an  M38  B-2  primer.  An  obturator  was  attached  to  the  projectile  during  firing  to  reduce 
combustion  gas  emissions.  Several  pretest  shots  verified  the  projectile  velocity  repeatability.  The  aim 
point  was  1.25  in  (31.75  mm)  below  the  fuze  top,  but  the  actual  impact  point  was  0.27  in  (6.8  mm)  below 
the  aim  point.  Based  on  the  results  of  the  simulation,  it  was  desired  to  have  a  radiographic  image  of  the 
impact  event  taken  250  ps  after  initial  impact.  The  actual  radiograph  was  taken  272  ps  after  impact, 
and  we  compared  this  radiograph  with  the  hydrocode  results  at  280  ps  (Figure  10).  At  this  instant,  the 
brass  core  of  the  fuze  is  bent  but  not  fractured,  and  the  aluminum  outerbody  appears  fractured  through 
the  cross  section.  The  projectile  deforms  by  slightly  bulging  beneath  the  point  of  impact,  and  the 
projectile's  tail  end  is  rotated  slightly  upward  at  this  instant  in  the  impact 

The  steel  projectile  penetrated  into  the  fuze  to  a  maximum  depth  of  approximately  13  mm  and 
formed  a  cireular  indentation  in  the  aluminum.  The  aluminum  subsequently  fractured,  and  the  fuze  top 
was  severed  (Figure  11).  The  fuze  top  fragment  also  contained  a  circular  impact  imprint  in  its  surface. 
Also  visible  in  Figure  11  is  the  bolt  that  holds  the  fuze  to  the  0.25-in  (6.35  mm)-thick  rolled- 
homogeneous  armor  (RHA)  base  and  the  failed  weld  line  in  the  RHA  base  plate.  The  weld  in  the  RHA 
base  plate  failed  during  impact,  and  the  entire  fuze  assembly  rotated  backwards  approximately  35 
degrees.  This  was  an  unanticipated  fixture  failure,  and  prevented  quantitative  comparisons  between  the 
experimental  and  numerical  fuze  deformation  results.  The  projectile  sustained  little  permanent 
deformation,  although  a  small  chip  was  visible  on  the  edge  beneath  the  impact  point. 

6.  COMPARISION  OF  NUMERICAL  AND  EXPERIMENTAL  RESULTS 

The  experimental  results  and  the  simulation  results  generally  agree  during  the  initial  stages  of  the 
impact  but  depart  significantly  at  later  stages  of  the  penetration.  The  reasons  for  the  observed  differences 
between  the  simulation  and  the  experiment  appear  in  the  discussion  following,  and  are  related  to:  1 )  the 
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inability  of  the  hydrocode  to  effectively  model  fracture,  and  2)  boundary  condition  differences  between 
the  experiment  and  those  assumed  in  the  numerical  simulation. 

The  most  obvious  difference  between  the  experimental  results  and  our  predictions  is  that  in  the 
experiment,  the  fuze  failed  by  fracture,  while  we  predicted  that  the  top  of  the  fuze  remained  intact  and 
connected  to  the  fuze  body  for  identical  impact  conditions.  Despite  the  fact  that  the  plastic  "hinge"  might 
suggest  a  failure  location,  the  inability  of  the  hydrocode  to  model  gross,  large-scale  fracture  in  structures 
is  a  primary  reason  why  there  is  a  disparity  between  numerical  impact  results  and  the  experiment. 
Instead,  the  hydrocode  simulates  fracture  failure  by  eliminating  an  element's  ability  to  sustain  shear  and 
tensile  normal  stress  when  the  level  of  damage  exceeds  a  critical  level  in  that  element.  To  address  this 
type  of  fracture  problem,  the  EPIC92  hydrocode  would  require  an  algorithm  that  permits  finite  element 
mesh  bifurcation  when  a  critical  energy  release  rate  is  attained  or  by  using  specialized  crack-tip 
singularity  finite  elements  (Ingraffea  1989). 

The  boundary  conditions  are  somewhat  different  between  the  experiment  and  the  numerical 
simulation,  despite  our  efforts  to  minimize  them.  Boundary  condition  differences  include:  1 )  modeling 
the  contact  surface  as  a  frictionless  interface,  whereas  in  reality,  contact  friction  can  influence  the 
deformation  in  low-speed  impact  problems,  2)  a  disparity  in  the  impact  point  and  impact  velocity 
between  the  experiment  and  the  simulation  (additional  hydrocode  simulations  could  have  been 
conducted  to  match  the  actual  impact  point  in  the  experiment;  however,  the  additional  time  and  expense 
needed  to  conduct  the  additional  simulations  would  not  have  added  further  insight  into  the  render-safe 
problem,  since  the  fundamental  fracture  physics  could  not  be  modeled  by  the  hydrocode),  and  3) 
interface  bonding  differences  between  the  brass  core  and  the  aluminum  outerbody.  Recall  that  in  the 
experiment,  the  brass  core  was  press-fit  into  the  aluminum  outerbody  and  that  given  enough  force,  slip 
could  occur  along  this  interface.  The  interface  was  modeled  so  as  to  prevent  slip  in  the  numerical 
simulation.  Furthermore,  failure  of  the  weld  line  in  the  experiment  permitted  the  fuze  body  to  rotate  at 
its  base,  whereas  the  base  of  the  fuze  was  rigidly  fixed  in  the  numerical  simulation. 
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7.  ANALYTICAL  APPROACH 


The  fuze-projectile  problem  addressed  earlier  is  reminiscent  of  a  cantilever  beam  that  is  subjected 
to  transverse  impact.  Although  a  complete  review  of  the  vast  literature  available  on  this  topic  is  beyond 
the  scope  of  this  report,  some  notable  works  relevant  to  the  current  problem  should  be  mentioned.  Much 
of  the  original  research  on  the  transverse  impact  of  elastic-plastic  beams  was  funded  by  organizations 
such  as  the  Office  of  Naval  Research  and  can  be  found,  for  example,  in  early  issues  of  publications  such 
as  the  ASME  Journal  of  Applied  Mechanics.  Both  Lee  and  Symonds  (1952)  and  Parkes  (1955) 
investigated  the  extent  of  applicability  of  using  the  methods  of  static  plasticity  solutions  to  dynamic 
problems.  Their  rigid-plastic  analyses  predict  the  development  of  a  moving  plastic  "hinge."  Interest¬ 
ingly,  a  plastic  hinge  develops  beneath  the  projectile  in  Event  2  (Figure  8),  a  deformation  feature  which 
is  common  in  problems  addressed  in  the  literature  involving  both  quasistatic  and  impulsively  loaded 
elastic-plastic  beams.  For  short  beams,  the  effect  of  shear  deformations  are  no  longer  negligible  and  their 
effects  can  be  estimated  by  construction  of  so-called  interaction  curves,  using  a  plane  stress  approach 
(Green  1954;  Drucker  1956)  or  using  a  variational  approach  (Hodge  1957).  Hodge's  (1957)  analysis 
shows  that  for  simply  supported  elastic-plastic  short  beams  in  which  beam  thickness  is  approximately 
equal  to  beam  length  (as  in  this  work),  the  maximum  load  that  the  beam  can  sustain  is  reduced  by 
approximately  20%  if  one  includes  the  effects  of  shear  in  the  analysis.  Some  of  the  discrepancies  between 
theory  and  experiment  discussed  in  Parkes'  (1955)  paper  were  reexamined  by  Ting  (1964, 1965),  whose 
analyses  rigorously  account  for  geometrical  effects  due  to  large  plastic  deformations  of  the  rigid-plastic 
beam.  Keer  and  Schonberg  (1986)  solve  an  elastic  cantilever  beam  problem  subjected  to  two  types  of 
indentation  displacement  fields  using  a  local  elasticity  solution  coupled  to  a  global  beam  theory  solution 
approach.  Both  the  local  deformation  effects  of  the  indentor  as  well  as  the  global  deformation  of  the 
cantilever  beam  are  modeled;  however,  the  solution  is  limited  to  elastic  media.  Finally,  Shu  et  al.  (1992) 
describe  the  mechanics  of  oblique  impact  of  impulsively  loaded  cantilever  beams  in  terms  of  axial  and 
flexural  components.  Their  model  includes  a  mass  of finite  size  at  the  beam  tip,  which  produces  a  double 
"hinge"  deformation  mechanism.  The  model  is  used  to  explain  why  the  observed  curvature  of  the  beam 
tip  is  less  than  that  predicted  by  Parkes’  (1955)  single-"hinge"  model.  Closed-form  three-dimensional 
analytical  solutions  for  transient  loading  of  composite  elastic-plastic  bodies  of  variable  cross  section  do 
not  appear  in  the  literature,  so  a  precise  comparison  with  the  hydrocode  results  is  not  possible.  A  closed- 
form  solution  to  the  elastic-plastic  short  beam  impact  problem  would  be  invaluable  to  the  explosive 
ordnance  disposal  community. 


Often,  the  analytical  solution  to  the  "real"  structural  mechanics  problem  is  intractable,  wherers 
useful  solutions  can  be  obtained  for  a  reduced  or  simplified  version  of  the  original  problem.  In  this  spirit, 
we  consider  the  problem  of  an  elastic  perfectly-plastic  cantilever  beam  of  circular  cross  section  subjected 
to  transverse  impact  as  a  model  of  the  fuze  impact  problem.  The  analysis  uses  elementary  beam  theory 
with  the  following  additional  simplifying  assumptions: 

1)  There  is  small  deformation  theory. 

2)  There  is  negligible  beam  inertia. 

3)  The  beam  material  is  homogeneous  (not  composite)  and  elastic  perfectly-plastic. 

4)  There  are  negligible  shear  stresses. 


The  remainder  of  this  section  is  devoted  to  developing  a  simple  analytical  cantilever  beam  model 
of  the  fuze  for  estimating  the  total  amount  of  strain  energy  absorbed  by  the  fuze  during  impact.  The 
analytical  model  results  are  then  compared  with  the  hydrocode  predictions. 


7.1  Deflections  of  a  Cantilever  Beam.  The  deflection  curve  of  a  beam  can  be  shown  to  be 
governed  by  the  following  second-order  ordinary  differential  equation  (Timoshenko  and  Gere  1972), 

-M=k  ,  (6) 

A?  El 
dx 

in  which  v  is  the  deflection,  M  is  the  bending  moment  applied  to  the  beam,  E  is  Young's  modulus,  I  is 
the  cross-sectional  moment  of  inertia  about  the  neutral  axis,  and  k  is  the  curvature  of  the  neutral  axis. 
The  quantity  El  is  known  as  the  flexural  rigidity  of  the  beam.  Two  successive  integrations  of  Equation 
(6),  with  boundary  conditions,  v(L)  =  v'(L)  =  0  (Figure  12a),  provide  the  solution,  v(x),  for  an  end-loaded 
cantilever  beam  as 


v(x)  = 
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in  which  L  is  the  beam  length,  P  is  the  magnitude  of  the  end  load,  and  the  bending  moment  is  given  by 
M  =  -Px.  The  maximum  deflection  occurs  at  the  cantilever  tip,  x  =  0,  and  is  given  by 
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12  Elastic-Plastic  Deflections  of  a  Cantilever  Beam.  Moment  equilibrium  of  stress  about  the 
neutral  axis  of  the  beam  yields  the  following  integral  expression  for  the  resultant  bending  moment  as 


M 


=  JI 


o  ydA  , 


(9) 


in  which  a  is  the  longitudinal  stress  acting  normal  to  an  element  of  cross-sectional  area,  dA,  and  y  is  the 
distance  from  the  neutral  axis  to  dA.  The  fundamental  assumption  in  slender  beam  theory  is  that  plane 
sections  remain  plane  and  normal  to  the  longitudinal  fibers  of  the  beam  so  that  the  longitudinal  strain 
e  =  K  y.  For  a  beam  made  of  an  elastic  (Hookean)  material  we  have,  a  =  E  e,  thus. 


o  =  K  E  y 


(10) 


It  is  clear  from  Equation  (10)  that  the  maximum  stress  occurs  at  the  outer  fibers  of  the  beam  (Figure  12b). 
For  an  elastic  perfectly-plastic  material  the  maximum  stress  is  the  yield  stress  (Figure  12c),  and  we  have 


a  =  K  E  e  , 
y 


(11) 


in  which  e  is  the  distance  from  the  neutral  axis  to  the  edge  of  the  elastic  core  (Figure  12b).  Substitution 
of  Equation  (10)  and  Equation  (11)  into  Equation  (9)  and  using  (d  A)  =  2  z  (dy)  in  Equation  (9)  for  a  beam 
of  circular  cross  section  and  diameter  h,  we  arrive  at 
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The  first  integral  on  the  right-hand  side  of  Equation  (12)  represents  the  contribution  to  the  bending 
moment  due  to  elastic  stresses,  while  the  second  integral  in  Equation  (12)  represents  the  contribution  to 
the  bending  moment  due  to  plastic  flow.  Yielding  begins  to  occur  in  the  circular  beam  when  M  = 
or  when 

lim  M  =  M  .  (14) 

e-»h/2  y 


which  gives  for  the  yield  moment,  M^,  for  the  circular  beam, 


M 

y 


,3 

o  h  7t 

y _ 
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(15) 


As  the  bending  moment  of  the  beam  increases,  the  region  of  purely  elastic  behavior  shrinks  until  the 
maximum  or  limiting  moment  in  the  beam  is  reached.  The  limit  is  known  as  the  plastic  moment ,  M^, 
and  is  obtained  as 

lim  M  =  M  (16) 

e  -40+  P 

For  the  beam  of  circular  cross  section  the  plastic  moment  is  given  by 


o  h 

%  =  V 

Timoshenko  and  Gere  (1972)  define  a  shape  factor,  f,  for  a  beam  of  given  cross-sectional  geometry  as 
the  ratio  of  the  plastic  moment  to  the  yield  moment,  or  f  =  M^/M^.  For  a  beam  of  circular  cross  section, 

f  =  16/371  =  1.70,  asobtainedbydividingEquation(17)byEquation(15).  In  the  linear-elastic  range, 
the  relation  between  bending  moment  and  curvature  can  be  nondimensionalized  (Timoshenko  and  Gere 
1972)  as  follows, 

=  -^  .  0<M<M  .  (18) 

M  y 

y  y 

Equation  (18)  is  a  linear  function  on  a  moment-curvature  diagram  (Figure  13).  However,  as  M 
increases  beyond  M^,  and  the  beam  begins  to  yield,  the  moment-curvature  relationship  becomes 
nonlinear  and  the  ordinate  value  approaches  an  asymptote  representing  the  plastic  moment,  M^,  which 
is  equivalent  to  the  shape  factor,  f,  as  shown  in  Figure  13.  In  the  nonlinear  range,  the  bending  moment 
in  Equation  (13)  is  nondimensionalized  by  dividing  by  Equation  (15)  to  obtain 
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in  which  k/k^  =  h/2e.  By  way  of  comparison,  it  can  be  shown  that  the  moment-curvature  relation  for 
a  rectangular  beam  takes  on  the  relatively  simple  form  (Timoshenko  and  Gere  1972), 

2 

^  =  M<M<M,  (20) 

My  2  2k2  P 

which  can  be  solved  for  the  curvature  ratio  in  terms  of  the  moment  ratio  as  follows. 


•The  presence  of  the  arcsine  function  in  Equation  (19)  makes  the  direct  algebraic  solution  for  k/k^ 
difficult.  However,  an  approximate  solution  may  be  obtained  by  plotting  Equation  (19)  forgiven  values 
of  K/Ky  and  curve  fitting  the  resulting  nonlinear  equation  using  a  Marquardt-Levenberg  (M-L)  regression 
algorithm  (Press  et  al.  1986).  Equation  (19)  is  plotted  (Figure  14)  together  with  the  following 
approximation  obtained  using  the  M-L  method. 

My  <  M  <  Mp  ,  (22) 


where  we  have  replaced  the  M-L  derived  numerical  value,  0.86810,  in  the  numerator  of  Equation  (22), 

with  V7/2  to  more  closely  resemble  the  form  of  the  numerator  in  Equation  (21).  Before  we  can  determine 
the  strain  energy  in  the  beam,  we  must  determine  the  tip  deflection,  8,  as  a  function  of  the  applied  load, 
P,  The  tip  deflection,  5,  of  the  elastic-plastic  beam  may  be  determined  by  use  of  the  second  curvature- 
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area  theorem  (Timoshenko  and  Gere  1972,  pp.  306-309)  as  follows, 


in  which  Xj  represents  the  limit  of  purely  elastic  behavior  in  the  beam.  The  tip  deflection  is  evaluated 
by  first  substituting  the  approximate  curvature  expression  (Equation  22)  into  the  integrand  of  the  second 
integral  in  Equation  (23),  and  using  the  identities,  =  M^/EI,  Xj=My/P,  and  My=PyL.  This  results  in  the 
following  approximate  expression  for  the  dimensionless  deflection  (plotted  in  Figure  15), 


3 

P  L 

in  which  8^  =  (from  Equation  (8))  is  the  maximum  elastic  deflection  at  the  yield  load  P^.  When 
P.=i6  ,  the  dimensionless  deflection  is  -^  =  2.55  . 

Py  3k  6^ 

7.3  Complementary  Energy  of  an  Elastic-Plastic  Cantilever  Beam.  The  dimensionless  comple¬ 
mentary  energy,  U‘,  in  a  beam  subjected  to  inelastic  bending  may  be  obtained  by  integrating 
Equation  (24)  with  respect  to  P/P^  as 


The  first  term  in  Equation  (25)  is  the  elastic  contribution  to  the  complementary  energy,  while  the  second 
term  is  the  elastic-plastic  contribution  to  the  complementary  energy.  Evaluation  of  these  integrals  gives  ■ 
the  total  dimensionless  complementary  energy  of  the  beam, 

u*  =  i  +  j?V3(f-pi)>^  +(Pi-l)  +V3VfT (3pi  -  2f -l)j/pi  ,  (26) 
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in  which  p,  =-£-  and  f  =  ii  .  At  the  maximum  bending  load,  p,  =f=i^  ,  and  u*  *1  +  1.005  =  1.505  in 

Py  3Jt  ^  *  3n  2 

Equation  (26).  Also, 

U*  +  U  =  -P-i-  ,  (27) 

P  8 
y  y 

in  which  U  is  the  dimensionless  strain  energy  of  the  beam.  Hence  in  Equation  (27),  U  =  (16/37t)(2.55) 
- 1 .505  =  2.82.  This  result  shows  that  the  the  total  dimensionless  strain  energy,  at  the  moment  the  plastic 
hinge  is  formed,  is  over  4.5  times  the  energy  absorbed  in  the  purely  elastic  bending  regime.  The  ratio 
of  plastic  bending  strain  energy  to  total  elastic  bending  strain  energy  is  thus  (2.82  -  8/37t)/(8/37c)  =  2.32. 

8.  COMPARISON  OF  NUMERICAL  AND  ANALYTICAL  RESULTS 


In  order  to  determine  the  strain  energy  in  the  elastic-plastic  beam,  we  must  first  calculate  realistic 
values  of  and  5^  using  the  material  properties  of  6061-T6  aluminum.  To  simplify  the  analysis,  the 
effect  of  the  brass  core  material  on  the  inelastic  bending  of  the  cantilever  beam  is  not  considered.  We 
assume  in  our  analysis  that  the  cantilever  beam  is  composed  entirely  of  aluminum  which  is  "stronger" 

than  the  brass  up  to  an  equivalent  plastic  strain  of  0.20  (Figure  3).  Therefore,  for  a  given  deflection,  our 

* 

analysis  would  overestimate  the  beam  bending  strain  energy  for  s  0.20,  and  underestimate  the  strain 

energy  for  e*’  >  0.20  .  However,  a  composite  beam  solution  could  be  obtained  by  generalizing 
Equation  ( 1 2)  to  account  for  contributions  to  the  moment  equilibrium  from  both  the  brass  and  aluminum 
materials.  The  load  at  the  onset  of  yield  is  a  strong  function  of  beam  thickness,  h,  and  is  given  by 


and  the  deflection  at  yield  is 


3 

M  o  h  jt 


3  2  2 

o  h  Ji  L  2o  L 

5  =  _y _ =  y 

y  96EI  3Eh 


(28) 


(29) 
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in  which  the  moment  of  inertia  i  =  .  With  a  =  3.25  kbars  (4.7  x  ICf  psi),  L  =  2.47  in  (62.7  mm) 

(measured  from  the  impact  point  to  the  base),  E  =  739  kbar  (1.07  x  10''  psi),  h  =  2.4  in  (61  mm). 
Equations  (28)  and  (29)  give,  =  25,  893  Ibf  (115  kN),  and  5^  =  0.0074  in  (0.188  mm). 


The  effect  of  beam  taper  can  be  analyzed  by  recognizing  that  the  fuze  geometry  consists  of  a 
nonprismatic  (tapered)  cone  in  which  the  diameter  can  be  approximated  with  the  expression  h  =  h(x)  = 
h^  +  (h^^-  h  Jx/L.  It  can  be  easily  shown  by  substitution  of  h(x)  into  Equation  (28)  that  the  load,  P^,  required 
to  cause  yield,  is  a  function  of  position  along  the  length  of  the  cone  of  diameter,  h^,  at  the  impact  location, 
and  of  diameter,  h^^,  at  the  built-in  end.  Therefore  Equation  (28)  can  be  rewritten  as 


P 


y 


32  X 


(30) 


with  similar  substitutions  required  in  the  expressions  for  beam  inertia,  I,  and  deflection,  5^.  In  a  prismatic 
cantilever  beam,  the  positions  of  maximum  stress,  and  bending  moment  are  coincident  and  occur 
at  the  built-in  end  with  diameter,  h^^.  However,  in  a  tapered  cantilever  beam,  we  find  that  the  maximum 
stress  is 


128  PL 


max 


2’  0  \{\  ■  K) 


at  X  = 


h  L 

a 


(31) 


in  which  the  position  of  is  determined  by  setting  dcfyidx  =  0  in  Equation  (30).  In  the  fuze  problem, 
hj^=  2.4  in  (61  mm)  and  h^  =  1.39  in  (35.3  mm),  so  that  occurs  at  x  =  0.689L,  measured  from  the 
impact  point  toward  the  built-in  end.  In  addition,  we  find  that  the  ratio  of  maximum  stress,  o  ,  to  the 

QUIX 

stress,  a^,  in  the  built-in  end  of  the  tapered  beam  to  be 

^max  = - ^  (32) 
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in  which,  for  the  specific  case  of  the  fuze  geometry,  we  find  that  is  5%  greater  at  x  =  0.689L  than 
at  X  =  L  (Oj^is  determined  by  letting  x=L  in  Equation  (30)).  For  the  specific  case  of  a  tapered  circular 
cantilever  beam  in  which  hj^=  2  h^,  Timoshenko  and  Gere  (1972)  find  a  19%  increase  in  relative  to 
which  occurs  at  x  =  0.5L.  Our  result  in  Equation  (32)  is  more  general,  however,  and  reduces  to  the 
Timoshenko  and  Gere  (1972)  result  when  hj^=  2  h^.  The  assumption  that  the  fuze  geometry  is  prismatic 
is  not  so  bad  since  there  is  only  a  5%  difference  between  and  in  a  beam  with  the  fuze  geometry 
considered  in  this  work. 

An  estimate  of  the  average  load,  P,  on  the  cantilever  beam  can  be  determined  from  the  length  of 
time  it  takes  the  projectile  to  decelerate  and  begin  to  rebound.  Inspection  of  the  fuze  centerline  nodal 
position  versus  time  plot  (Figure  16)  indicates  that  the  projectile's  velocity  reversed  at  about  360  |is  and 
that  the  amount  of  elastic  rebound  beneath  the  projectile  is  about  0.0298  in  (0.747  mm).  The  change  in 
the  projectile's  momentum  with  respect  to  time  is  equal  to  the  average  impact  force  P  =  37,  094  Ibf 
(165  kN),  which  is  determined  using  an  initial  projectile  velocity  of  198  m/s  and  mass  of  0.3  kg.  The 
beam  theory  solution  provided  the  value  of  P^  =  25,  893  Ibf  (115  kN)  from  Equation  (28),  so  that 
P/P^  =  1 .43.  Since  P^  is  sensitive  to  the  value  chosen  for  h  (see  Equation  (28)),  if  a  slightly  smaller  value 
of  h  were  used  in  the  calculation  (rather  than  the  maximum  value  used  of  h^= h  =  2.4  (61  mm)),  P^  would 
•dramatically  decrease,  and  P/P^,  would  become  f  =  M^/M^  =  16/371  =  1.7.  At  this  limit  the  beam  is  in  a 
state  of  unrestricted  plastic  flow.  A  state  of  unrestricted  plastic  flow  occurs  in  the  fuze  very  early  in  the 
deformation  history,  as  seen  in  the  well-developed  pattern  of  plastic  strain  through  the  fuze  cross  section 
at  120  ps  (Figure  17).  The  total  strain  energy  absorbed  by  the  fuze  during  impact  can  be  estimated  by 
assuming  that  P/P^  =  M^/M^  =  1 6/37t  =  1 .7  and  5^  =  0.0074  in  (0. 1 88  mm)  from  Equation  (29).  The  actual 
deflection  beneath  the  load  point  from  the  numerical  simulation  at  360  ps  is  5  =  0.657  in  (16.7  mm), 
which  results  in  5/6^  =  89.  Hence,  the  total  dimensionless  strain  energy  absorbed  by  the  fuze  due  to 
elastic-plastic  bending  is  estimated  with  the  aid  of  Figure  15  to  be 

=  -L(A-2.55)  +  U=  -1^89-2.55)  +  2.8=  149.5  .  (33) 

P  0  3jc 

y  y 

This  energy  is  dimensionalized  through  multiplication  by  PyS^,  which  results  in  a  total  strain  energy  of 

2,387  ft-lbf  or  3.24  kJ,  attributed  to  elastic-plastic  bending  of  the  fuze.  The  0.3-kg  projectile  can  impart 
5.88  kJ  of  kinetic  energy  during  impact.  Therefore,  the  plastic  bending  energy  is  55%  of  the  initial 
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projectile  kinetic  energy.  The  fuze  also  suffered  a  considerable  amount  of  localized  deformation  beneath 
the  projectile  (Figure  17).  The  energy  absorbed  by  this  localized  deformation  can  be  approximated  by 
multiplying  the  average  projectile  indentation  depth,  5^  (estimated  from  the  experiment),  with  the 
average  impact  force,?,  resulting  in  =(37,094)(0.45)inlbf=  1.89kJ.  This  localized  plastic  energy 
is  32%  of  the  initial  projectile  kinetic  energy.  Hence,  the  energy  absorbed  in  fuze  bending  is  about  1 .7 
times  greater  than  the  energy  dissipated  by  localized  deformation  beneath  the  projectile.  In  the 
Event  2  simulation,  the  EPIC92  code  estimates  that  the  total  plastic  work  absorbed  by  the  fuze  is  91% 
of  the  initial  projectile  kinetic  energy,  which  is  in  close  agreement  to  our  analysis,  which  predicts  that 
87%  of  the  initial  projectile  kinetic  energy  is  dissipated  in  plastic  work.  The  EPIC92  code  also  predicts 
that  the  rotational  kinetic  energy  associated  with  the  rebounding  projectile  is  only  1%  of  the  initial 
projectile  kinetic  energy. 

The  analytical  model  of  an  elastic-plastic  beam  subjected  to  an  end  load  provides  insight,  at  least 
to  the  point  of  fracture,  into  the  render-safe  mechanics  of  fuzes.  It  is  estimated  that  55%  of  the  initial 
projectile  kinetic  energy  is  dissipated  in  elastic-plastic  fuze  bending  and  32%  of  the  energy  is  dissipated 
in  localized  deformation  beneath  the  projectile.  These  results  are  corroborated  by  the  EPIC92  numerical 
results,  which  predict  that  91  %  of  the  projectile  kinetic  energy  is  dissipated  into  target  plastic  work;  the 
remaining  energy  is  primarily  dissipated  as  plastic  work  in  deforming  the  projectile. 

9.  CONCLUSIONS 

1)  The  results  of  numerical  hydrocode  simulations  of  oblique  impact  of  a  cylindrical  steel 
projectile  onto  an  aluminum-brass  composite  fuze  simulant  indicate  that  projectile  hardness  is  an 
important  parameter  that  governs  the  deformation  history  of  the  fuze;  a  "soft"  steel  projectile  impacted 
the  fuze  and  created  a  gouge,  elliptical  in  form,  as  it  ricocheted  off  the  fuze,  whereas  a  "hard"  steel 
projectile  penetrated  to  a  greater  depth  and  bounced  obliquely  off  the  fuze. 

2)  Reasons  for  the  observed  disparity  between  the  experiment  and  the  hydrocode  simulation  of 
a  materially  and  geometrically  identical  fuze  impact  problem  include:  1)  the  inability  of  the  hydrocode 
to  effectively  model  global  fracture,  and  2)  boundary  condition  differences  between  the  experiment  and 
the  numerical  simulation;  in  the  experiment  a  weld  line  failed,  which  provided  an  additional  rotational 
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degree  of  freedom,  whereas  in  the  numerical  simulation  the  fuze  maintained  a  rigid  base.  Nonetheless, 
the  plastic  "hinge"  predicted  in  the  simulation  might  suggest  an  indication  of  where  the  failure  might  be 
expected  to  occur. 

3)  The  projectile's  sense  of  rotation  during  impact  is  related  to  projectile  hardness  and  impact 
point  as  the  numerical  simulations  indicate  that  the  projectile  1)  rotates  clockwise  in  impact  Event  1, 
involving  the  "soft"  projectile,  because  of  slip  between  the  projectile  and  fuze,  2)  rotates  counterclock¬ 
wise  in  impact  Event  2,  involving  the  "hard"  projectile,  since  the  harder  projectile  penetrates  into  the  fuze, 
thereby  precluding  slip,  and  3)  remains  relatively  level  in  impact  Event  3  as  the  obliquity  of  impact  is 
less  than  in  Events  1  and  2. 

4)  Positioning  the  gun  close  to  the  fuze  during  firing  provides  a  lateral  constraint  to  the  projectile 
during  impact.  This  practice  would  help  prevent  projectile  rotation  during  flight  and  impact  and  would 
also  reduce  the  variability  in  impact  response  that  occurs  as  a  result  of  projectile  impact  obliquity,  impact 
location,  and  projectile  hardness  effects. 

5)  An  analytical  model  of  a  cantilevered  elastic-plastic  beam  subjected  to  an  end  load  provides 
insight,  at  least  to  the  point  of  fracture,  into  the  render-safe  mechanics  of  fuzes.  For  the  case  studied, 
it  is  estimated  that  55%  of  the  initial  projectile  kinetic  energy  is  dissipated  in  elastic-plastic  fuze  bending 
and  32%  of  the  energy  is  dissipated  in  localized  deformation  beneath  the  projectile.  These  results  are 
corroborated  by  the  EPIC92  numerical  results,  which  predict  that  91%  of  the  projectile  kinetic  energy 
is  dissipated  into  target  plastic  work;  the  remaining  energy  is  primarily  dissipated  as  plastic  work  in 
deforming  the  projectile. 
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LIST  OF  SYMBOLS 


D  EPIC  damage  parameter 

dA  beam  area  increment 

Djjlj  elastic  modulus  tensor 

E  Young’s  modulus 

e  limit  of  elastic  behavior 

h  beam  diameter 

tapered  beam  diameter  at  tip 
tapered  beam  diameter  at  built-in  end 
I  beam  moment  of  inertia 

L  beam  length 

M  beam  moment 

M^.  beam  yield  moment 

M  beam  plastic  moment 

P load 

P^,  load  at  initial  yield 

T’  homologous  temperature 

U  dimensionless  strain  energy 

U’  dimensionless  complementary  energy 

total  dimensionless  strain  energy 

V  elastic  beam  deflection 

V  maximum  elastic  tip  deflection 

X  beam  axial  coordinate 

* 

Xj  limit  of  elastic  behavior 

y  beam  coordinate 

5  beam  tip  deflection 

5^,  tip  deflection  at  initial  yield 

da ,  incremental  stress  tensor 

de..  incremental  strain  tensor 

ij 

de^  incremental  plastic  strain  tensor 

e  strain  rate 

e  o  normalization  strain  rate 

.  * 

e  dimensionless  plastic  strain  rate 

fracture  strain 
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equivalent  plastic  strain 

Ae^  incremental  equivalent  plastic  strain 

K  beam  curvature 

K^,  beam  curvature  at  initial  yield 

a  beam  stress 

o  equivalent  stress 

* 

a  pressure -stress  ratio 

o  mean  normal  stress 

a  maximum  beam  stress 

rrux 

maximum  bending  stress 
o' .  deviatoric  stress 


38 


NO.  OF 

COPIES  ORGANIZATION 

2  ADMINISTRATOR  DTIC 

ATTN  DTIC  DDA 
CAMERON  STATION 
ALEXANDRIA  VA  22304-6145 

1  CDR  USAMC 

ATTN  AMCAM 
5001  EISENHOWER  AVE 
ALEXANDRIA  VA  22333-0001 

1  DIR  USARL 

ATTN  AMSRL  OP  SD  TA 
2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1145 

3  DIR  USARL 

ATTN  AMSRL  OP  SD  TL 
2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1145 

1  DIR  USARL 

ATTN  AMSRL  OP  SD  TP 
2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1145 

2  CDR  US  ARMY  ARDEC 
ATTN  SMCAR  TDC 
PCTNY  ARSNL  NJ  07806-5000 

1  DIR  BENET  LABS 

ATTN  SMCAR  CCB  TL 
WATERVLIETNY  12189-4050 

1  DIR  USA  ADVANCED  SYSTEMS 

R&A  OFC 

ATTN  AMSAT  R  NR  MS  219  1 
AMES  RESEARCH  CENTER 
MOFFETT  FLD  CA  94035-1000 

1  CDR  US  ARMY  MICOM 

ATTN  AMSMl  RD  CS  R  DOC 
RDSTN  ARSNL  AL  35898-5010 

1  CDR  US  ARMY  TACOM 

ATTN  AMSTA  JSK 
ARMOR  ENG  BR 
WARREN  MI  48397-5000 


NO.  OF 

COPIES  ORGANIZATION 
1  DIR 

USA  TRADOC  ANALYSIS  CMD 
ATTN  ATRC  WSR 
WSMR  NM  88002-5502 

1  CMDT 

US  ARMY  INFANTRY  SCHOOL 
ATTN  ATSH  CD  SECURITY  MGR 
FT  BENNING  GA  31905-5660 


ABERDEEN  PROVING  GROUND 

2  DIR  USAMSAA 

ATTN  AMXSY  D 

AMXSY  MP  H  COHEN 

1  CDR  USATECOM 

ATTN  AMSTE  TC 

1  DIR  USAERDEC 

ATTN  SCBRD  RT 

1  CDR  USACBDCOM 

ATTN  AMSCB  CH 

1  DIR  USARL 

ATTN  AMSRL  SL  1 

5  DIR  USARL 

ATTN  AMSRL  OP  AP  L 


39 


NO.  OF 

COPIES  ORGANIZATION 

1  HQDA  SARD  TT  DR  F  MILTON 
WASH  DC  20310-0103 

1  HQDA  SARD  TT  MR  J  APPEL 
WASH  DC  20310-0103 

1  HQDA  SARD  TT  MS  C  NASH 
WASH  DC  20310-0103 

1  HQDA  SARD  TR  DR  R  CHAIT 

WASH  DC  20310-0103 

1  HQDA  SARD  TR  MS  K  KOMINOS 

WASH  DC  201310-0103 

1  DIR  USARL 

ATTN  AMSRL  CP  CA  D  SNIDER 
2800  POWDER  MILL  ROAD 
ADELPHl  MD  20783-1145 

8  DIR  USARL 

ATTN  AMSRL  MA  P 
L  JOHNSON 
B  HALPIN 
T  CHOU 
C  WHITE 
J  MCLAUGHLIN 
AMSRL  MA  PA 
D  GRANVILLE 
W'  HASKELL 

AMSRL  MA  MA  G  HAGNAUER 
WATERTOWN  MA  02172-0001 

4  CDR  US  ARMY  ARDEC 
ATTN  SMCAR  FSE 
T  GORA 

E  ANDRICOPOULOS 
B  KNTJTELSKY 
A  GRAF 

PCTNY  ARSNTL  NJ  07806-5000 

3  CDR  US  ARMY  ARDEC 

ATTN  SMCAR  TD 
R  PRICE 
V  LINDER 
T  DAVIDSON 

PCTNW  ARSNT  NJ  07806-5000 

1  CDR  US  ARMY  ARDEC 

ATTN  F  MCLAUGHLIN 
PCTNY  ARSNL  NJ  07806 


NO.  OF 

COPIES  ORGANIZATION 

5  CDR  US  ARMY  ARDEC 
ATTN  SMCAR  CCH  T 
S  MUSALLI 
P  CHRISTIAN 
K  FEHSAL 
N  KRASNOW 
R  CARR 

PCTNY  ARSNL  NJ  07806-5000 

1  CDR  US  ARMY  ARDEC 

ATTN  SMCAR  CCH  V  E  FENNELL 
PCTNY  ARSNL  NJ  07806-5000 

1  CDR  US  ARMY  ARDEC 

ATTN  SMCAR  CCH  J  DELOREN^O 
PCTNY  ARSNL  NJ  07806-5000 

2  CDR  US  ARMY  ARDEC 
ATTN  SMCAR  CC 

J  HEDDERICH 

COL  SINCLAIR 

PCTN^  ARSNL  NJ  07806-5000 

1  CDR  US  ARMY  ARDEC 
ATTN  SMCAR  CCH  P  J  LUTZ 
PCTNT  ARSNL  NJ  07806-5000 

2  CDR  US  ARMY  ARDEC 
ATTN  SMCAR  FSA  M 
D  DEMELLA 

F  DIORIO 

PCTNY  ARSNL  NJ  07806-5000 

1  CDR  US  ARMY  ARDEC 

ATTN  SMCAR  FSA  C  SPINELLI 
PCTNY  ARSNL  NJ  07806-5000 

11  DIR  BENET  LABS 

ATTN  SMCAR  CCB 
C  KITCHENS 
J  KEANE 
J  SANTINI 
J  VASILAKIS 
G  FRIAR 
R  nSCELLA 
V  MONTVORI 
J  WRZOCHALSKI 
J  BATTAGLIA 
R  HASENBEIN 
SMCAR  CCB  R  S  SOPOK 
WATERVLIETNY  12189 


40 


NO.  OF 

COPIES  ORGANIZATION 

1  CDR  WATERVLIET  ARSENAL 

ATTN  SMCWV  QAE  Q  C  HOWD 
BLDG  44 

WATERVLIET  NY  12189-4050 

1  CDR  WATERVLIET  ARSENAL 

ATTN  SMCWV  SPM 
T  MCCLOSKEY 
BLDG  25  3 

WATERVLIET  NY  12189-4050 
1  CDR 

WATERVLIET  ARSENAL 
ATTN  SMCWV  QA  QS  K  INSCO 
WATERVLIET  NY  121894050 

1  CDR  PRODUCTION  BASE 

MODERNIZATION  ACTIVITY 
ATTN  AMSMC  PBM  K 
PCTNY  ARSNL  NJ  07806-5000 

1  CDR  USABRDC 

ATTN  STRBE  JBC  C  KOMINOS 
FT  BELVOIR  VA  22060-5606 

1  DIR  USACRREL 

ATTN  P  DUTTA 
72  LYME  ROAD 
HANOVER  NH  03755 

1  DIR  USARL 

ATTN  AMSRL  WT  L  D  WOODBURY 
2800  POWDER  MILL  ROAD 
ADELPHI  MD  20783-1145 

4  CDR  US  ARMY  MICOM 

ATTN  AMS  Ml  RD  W  MCCORKLE 
AMS  MI  RD  ST  P  DOYLE 
AMS  MI  RD  ST  CN  T  VANDIVER 
AMSMI  RD  ST  WF  MIKE  COLE 
RDSTN  ARSNL  AL  35898 

2  US  ARMY  RESEARCH  OFFICE 
MATH  &  COMPUTER  SCIENCES  DIV 
ATTN  ANDREW  CROWSON 

J  CHANDRA 

PO  BOX  12211 

RSCH  TRI  PK  NC  27709-2211 


NO.  OF 

COPIES  ORGANIZATION 

2  US  ARMY  RESEARCH  OFFICE 
ENGINEERING  SCIENCES  DIV 
ATTN  G  ANDERSON 
R  SINGLETON 
PO  BOX  12211 

RSCH  TRI  PK  NC  27709-2211 

2  PM  SADARM 

PCTNY  ARSNL  NJ  07806-5000 

2  PM  TMAS 

ATTN  SFAE  AR  TMA 
COL  BREGARD 
C  KIMKER 

PCTNY  ARSNL  NJ  07806-5000 

3  PM  TMAS 

ATTN  SFAE  AR  TMA  MD 
H  YUEN 
J  MCGREEN 
R  KOWALSKI 

PCTNY  ARSNL  NJ  07806-5000 

2  PM  TMAS 

ATTN  SFAE  AR  TMA  MS 
R  JOINSON 
D  GUZIEWICZ 

PCTNY  ARSNL  NJ  07806-5000 

1  PM  TMAS 

ATTN  SFAE  AR  TMA  MP  W  LANG 
PCTNY  ARSNL  NJ  07806-5000 

2  PEO  ARMAMENTS 
ATTN  SFAE  AR  PM 
D  ADAMS 

T  MCWILLIAMS 

PCTNY  ARSNL  NJ  07806-5000 

1  PEO  FIELD  ARTILLERY  SYSTEMS 

ATTN  SFAE  FAS  PM  H  GOLDMAN 
PCTNY  ARSNL  NJ  07806-5000 

4  PM  AFAS 

ATTN  LTC  D  ELLIS 
G  DELCOCO 
J  SHIELDS 
B  MACHAK 

PCTNY  ARSNL  NJ  07806-5000 


41 


NO.  OF 

COPIES  ORGANIZATION 


NO.  OF 

COPIES  ORGANIZATION 


2  CDR  WRIGHT  PATTERSON  AFB 
ATTN  WL  MLBM  S  DONALDSON 
WL  nv  A  MAYER 
DAYTON  OH  45433-7750 

2  NASA  LANGLEY  RSCH  CTR 
MAIL  STOP  266 
ATTN  AMSRL  VS  W  ELBER 
AMSRL  VS  S  F  BARTLETT  JR 
HAMPTON  VA  23681-0001 

2  NAVAL  SURFACE  WARFARE  CTR 
CODE  G33 

DAHLGREN  VA  22448 

1  OFFICE  OF  NAVAL  RESEARCH 

MECH  DIV  CODE  1132SM 
ATTN  YAP  A  RAJAPAKSE 
ARLINGTON  VA  22217 

1  NAVAL  ORDNANCE  STATION 
ADV  SYS  TECHNOLOGY  BR 
ATTN  D  HOLMES  CODE  2011 
LOUISVILLE  KY  40214-5245 

1  DAVID  TAYLOR  RSCH  CTR 
SHIP  STRCTS  &  PRTCTN  DEPT 
ATTN  J  CORRADO  CODE  1702 
BETHESDA  MD  20084 

2  DAVID  TAYLOR  RSCH  CTR 
ATTN  R  ROCKWLLL 

W  PHYILLAIER 
BETHESDA  MD  20054-5000 

5  DIR  LLNL 

ATTN  R  CHRISTENSEN 
S  DETERESA 
W  FENG 
F  MAGNESS 
M  HNGER 
PO  BOX  808 
LIVERMORE  CA  94550 

1  DIR  LANL 

ATTN  D  RABERN 
MEE  13  MAIL  STOP  J  576 
PO  BOX  1633 
LOS  ALAMOS  NM  87545 


1  OAK  RIDGE  NATIONAL  LAB 
ATTN  R  M  DAVIS 

PO  BOX  2008 

OAK  RIDGE  TN  37831-6195 

2  BATTELLE  PNL 
ATTN  M  SMITH 
M  C  C  HAMPTON 
PO  BOX  999 
RICHLAND  WA  99352 

6  DIR  SNL 

ATTN  C  ROBINSON 
G  BENEDETTI 
W  KAWAHARA 
K  PERANO 
D  DAWSON 
PNIELAN 
PO  BOX  969 

LIVERMORE  CA  94550-0096 

1  PENTI  STATE  UNIV 

ATTN  RICHARD  MCNITT 
227  HAMMOND  BLDG 
UNIVERSITY  PARK  PA  16802 

1  SOUTHWEST  RSCH  INSTITUTE 

ATTN  C  ANDERSON 
6220  CULEBRA  ROAD 
SAN  ANTIONIO  TX  78284 

1  UCLA 

MANE  DEPT  ENGRG  IV 

ATTN  H  THOMAS  HAHN 
LOS  ANGELES  CA  90024-1597 

2  UNIV  OF  DAYTON  RSCH  INST 
ATTN  RAN  Y  KIM 

AJIT  K  ROY 

300  COLLEGE  PARK  AVE 
DAYTON  OH  45469-0168 

2  ■  UNIV  OF  DELAWARE 

CTR  FOR  COMPOSITE  MTRLS 
ATTN  J  GILLESPE 
M  SANT ARE 

201  SPENCER  LABORATORY 
NEWARK  DE  19716 


42 


NO.  OF 

COPIES  ORGANIZATION 

1  UNIV  OF  TEXAS  AT  AUSTIN 
CTR  FOR  ELECTROMECHANICS 
ATTN  J  PRICE 
10100  BURNET  ROAD 
AUSTIN  TX  787584497 

1  AAI  CORPORATION 
ATTN  TECH  LIBRARY 
PO  BOX  126 

HUNT  VALLEY  MD  21030-0126 

1  ARMTEC  DEFENSE  PRODUCTS 
ATTN  STEVE  DYER 
PO  BOX  848 
COACHELLA  CA  92236 

3  ALLIANT  TECHSYSTEMS  INC 
ATTN  J  BODE 
C  CANDLAND 
K  WARD 

5901  LINCOLN  DR 
MINNEAPOLIS  MN  55346-1674 

1  ALLIANT  TECHSYSTEMS  INC 

ATTN  TIM  HOLMQUIST 
600  SECOND  STREET  NE 
HOPKINS  MN  55343 

1  BALLISTIC  IMPACT  DYNAMICS 

ATTN  RODNEY  RECHT 
3650  S  CHEROKEE  2 
ENGLEWOOD  CO  80110 

1  CALIFORNIA  RSRCH  &  TECH 

ATTN  DENNIS  ORPHAL 
5117  JOHNSON  DRIVE 
PLEASANTON  CA  94566 

1  CHAMBERLAIN  MFG  CORP 

R&D  DIVISION 
ATTN  M  TOWNSEND 
PO  BOX  2545 
WATERLOO  lA  50704 

1  COMP  MECHANICS  ASSOC 

ATTN  JONAS  A  ZUKAS 
PO  BOX  11314 
BALTIMORE  MD  21239-0314 


NO.  OF 

COPIES  ORGANIZATION 

1  CUSTOM  ANALYTICAL  ENGRNG 
SYSTEMS  INC 
ATTN  A  ALEXANDER 
STAR  ROUTE  BOX  4A 
FLINTSTONEMD  21530 

1  GDLS  DIVISION 
ATTN  D  BARTLE 
PO  BOX  1901 
WARREN  MI  48090 

3  HERCULES  INCORPORATED 
ATTN  G  KUEBELER 
J  VERMEYCHUK 
B  MANDERVILLE  JR 
HERCULES  PLAZA 
WILMINGTON  DE  19894 

1  HEXCEL 

ATTN  M  SHELENDICH 
PO  BOX  2312 
DUBLIN  CA  94568-0705 

1  lAP  RESEARCH  INC 
ATTN  A  CHALLITA 
2763  CULVER  AVENUE 
DAYTON  OHIO  45429 

5  INSTITUTE  FOR  ADVANCED 
TECHNOLOGY 
ATTN  T  KIEHNE 
HFAIR 
P  SULLIVAN 
S  BLESS 

R  SUBRAMANIAN 
4030  2  W  BRAKER  LANE 
AUSTIN  TX  78759 

1  INTERFEROMETRICS  INC 
ATTN  R  LARRIVA  VP 
8150  LEESBURG  PIKE 
VIENNA  VA  22100 

3  KAMAN  SCIENCES  CORP 
ATTN  D  ELDER 
T  HAYDEN 
N  ARI 

PO  BOX  7463 

COLORADO  SPRINGS  CO  80933 


43 


NO.  OF 
COPIES 

1 


1 


2 


1 


2 


2 


1 


2 


ORGANIZATION 


NO.  OF 

COPIES  ORGANIZATION 


D  R  KENNEDY  &  ASSOC  INC 
ATTN  D  KENNEDY 
PO  BOX  4003 

MOUNTAIN  VIEW  CA  94040 

LIVERMORE  SOFTWARE  TECH 
CORP 

ATTN  J  O  HALLQUIST 
2876  WAVERLY  WAY 
LIVERMORE  CA  94550 

LORAL  VOUGHT  SYS  CORP 
ATTN  G  JACKSON 
K  COOK 

1701  3\TST  MARSHALL  DR 
GRANT)  PRAIRIE  TX  75051 

LORAL  VOUGHT  SYS  CORP 
ATTN  KEN  HAVENS  MS  EM  36 
PO  BOX  650003 
DALLAS  TX  75265-0003 

MARTIN  MARIETTA  CORP 
ATTN  P  DEV'AR 
L  SPONAR 

230  EAST  GODDARD  BLVD 
KING  OF  PRUSSIA  PA  19406 

OLIN  CORPORATION 
FLINCHBAUGH  DIV 
ATTN  E  STEINTR 
B  STEWART 
PO  BOX  127 
RED  LION  PA  17356 

OLIN  CORPORATION 
ATTN  L  WHITMORE 
10101  9TH  ST  NORTH 
ST  PETERSBURG  FL  33702 

UNIV  OF  MINNESOTA  AHPCRC 
ATTN  G  SELL 
D  AUSTIN 

1100  WASHINGTON  AVE  S 
MINNEAPOLIS  MN  55415 


ABERDEEN  PROVING  GROUND 


109  DIR  USARL 

ATTN  AMSRL  Cl  C  MERMAGEN  394 
AMSRL  Cl  C  W  STUREK  1121 
AMSRL  Cl  CB  R  KASTE  394 
AMSRL  Cl  S  A  MARK  309 
AMSRL  SL  B  P  DEITZ  328 
AMSRL  SL  BA  J  WALBERT  1065 
AMSRL  SL  BL  D  BELY  328 
AMSRL  SL  I  D  HASKHJL  1065 
AMSRL  WT  P  A  HORST  390A 
AMSRL  WT  PA 
T  MINOR  390 
C  LEVERITT  390 
D  KOOKER  390A 
AMSRL  WT  PB 
E  SCHMIDT  120 
P  PLOSTINS  120 
AMSRL  WT  PC  R  FIFER  390A 
AMSRL  WT  PD 
B  BURNS  390 
W  DRYSDALE  390 
K  BANNISTER  390 
T  BOGETTI  390 
J  BENDER  390 
R  MURRAY  390 
R  KIRKENDALL  390 
TERLINE  390 
D  HOPKINS  390 
S  WILKERSON  390 
C  MCCALL  390 
D  HENRY  390 
R  KASTE  390 
L  BURTON  390 
J  TZENG  390 
R  LIEB  390 
G  GAZONAS  390  5  CP 
M  LEADORE  390 
AMSRL  WT  PD  ALC 
A  ABRAHAMIAN 
K  BARNES 
M  BERMAN 
H  DAVISON 
AFRYDMAN 
TLI 

W  MCINTOSH 
E  SZYMANSKI 


44 


NO.  OF 
COPIES 


ORGANIZATION 


NO.  OF 

COPIES  ORGANIZATION 


AMSRL  WT  T 

W  MORRISON  309 
T  WRIGHT  309 
AMSRL  WT  TA 
W  GILLICH  390 
W  BRUCHEY  390 
J  DEHN  390 
AMSRL  WT  TB 
K  BENJAMIN  309 
T  DORSEY  309 
R  FREY  309 
F  GREGORY  309 
W  HILLSTROM  309 
W  LAWRENCE  309 
O  LYMAN  1185 
J  STARKENBERG  309 
L  VANDE  KIEFT  309 
J  WATSON  309 
V  BOYLE  309  5  CP 
S  STEGALL  1185  5  CP 
AMSRL  WT  TC 

W  DE  ROSSET  309 
K  KIMSEY  309 
R  COATES  309 
F  GRACE  309 
M  LAMPSON  309 
D  SCHEFFLER  309 
B  SORENSEN  309 
R  SUMMERS  309 
E  WALKER  309 
W  WALTERS  309 
AMSRL  WT  TD 
D  DIETRICH  309 
G  RANDERS  PEHRSON  309 
J  HUFFINGTON  309 
A  DAS  GUPTA  309 
J  SANTIAGO  309 
K  FRANK 
J  HARRISON 
M  SCHEIDLER 
S  SEGLETES  5  CP 
J  WALTER 
AMSRL  WT  NC 
R  LOTTERO  309 
R  PEARSON  309 
S  SCHRAML  309 
AMSRL  WT  W  C  MURPHY  120 
AMSRL  WT  WA 
H  ROGERS  394 
B  MOORE  394 
A  BARAN  394 


AMSRL  WTWB 
F  BRANDON  120 
W  D’AMICO  120 
AMSRL  WT  WC  J  ROCCHIO  120 
AMSRL  WTWD 
A  ZIELINSKI  120 
J  POWELL  120 
A  PRAKASH  120 
AMSRL  WT  WE 
J  TEMPERLEY  120 
J  THOMAS  394 


45 


Intentionally  left  blank. 


46 


USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 

This  Laboratory  undertakes  a  continuing  effort  to  improve  the  quality  of  the  reports  it  publishes.  Your 
comments/answers  to  the  items/questions  below  will  aid  us  in  our  efforts. 

1.  ARL  Report  Number  ARl-TR-665 _ Date  of  Report _ January  1995 _ 

2.  Date  Report  Received _ 

3.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or  other  area  of  interest  for 

which  the  report  will  be  used.) _ _ 


4.  Specifically,  how  is  the  rejxirt  being  used?  (Information  source,  design  data,  procedure,  source  of 
ideas,  etc.) _ 


5.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far  as  man-hours  or  dollars  saved, 
operating  costs  avoided,  or  efficiencies  achieved,  etc?  If  so,  please  elaborate. _ 


6.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future  reports?  (Indicate 
changes  to  organization,  technical  content,  format,  etc.) _ 


Organization 

CURRENT  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 

City,  State,  Zip  Code 

7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the  Current  or  Correct  address 
above  and  the  Old  or  Incorrect  address  below. 

Organization 

OLD  Name 

ADDRESS  _ _ _ 

Street  or  P.O.  Box  No. 

City,  State,  Zip  Code 
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